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Preface 


» 

This  report  presents  the  results  of  my  investigation 

into  three  methods  for  obtaining  the  optimal  control  re- 
quirements of  satellites  with  flexible  appendages.  A new 
method  using  integral  coordinate  techniques  was  introduced 
in  this  thesis  and  was  found  to  provide  quite  acceptable 
optimal  control  results. 

This  project  required  that  I become  fairly  involved 
with  the  concepts  of  modern  control  theory,  stability  anal- 
ysis, and  dynamics  of  flexible  vehicles.  In  this  regard, 
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areas.  I am  especially  grateful  to  Dr.  Robert  A.  Calico, 
my  thesis  advisor,  for  suggesting  the  original  idea  for 
this  thesis  and  for  the  guidance  he  provided  throughout 
this  project. 

Finally,  I would  like  to  thank  my  wife,  Maria,  and 
i our  daughters,  Maria  Angela,  Kristina,  and  Diana  for  their 

patience  and  understanding  during  the  past  eighteen  months. 
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Abstract 


•^^This  thesis  investigated  three  methods  for  obtaining 
optimal  control  requirements  for  satellites  having  flexible 
appendages.  A discrete  method,  modal  method,  and  a new 
integral  coordinate  method  for  obtaining  control  require- 
ments were  examined.  Mathematical  formulation  was  performed 
on  a satellite  configuration  consisting  of  a symmetrical 
rigid  body  with  two  flexible  antennas  extending  in  opposite 
directions  along  the  spin  axis.  System  equations  of  motion 
were  derived  using  Hamilton's  equations.  Modern  optimal 
control  theory,  involving  the  minimization  of  a quadratic 
cost  functional  and  the  numerical  solution  to  the  steady 
state  matrix  Riccati  equation,  was  applied  to  the  system. 

An  example  problem  was  presented  and  numerically  solved. 

The  resulting  controlled  states  and  control  requirements 
for  each  method  were  compared  and  a discussion  of  the 
advantages  and  disadvantages  of  the  different  techniques 
was  presented.  J^The  results  of  this  thesis  indicate  that 
the  integral  coordinate  technique  provides  a valid  and  use- 
ful means  of  obtaining  realistic  estimates  of  control  re- 
quirements for  the  class  of  satellites  under  consideration. 

r- 

The  report  goes  on  to  recommend  further  investigation  of 
the  integral  coordinate  method  for  other  spacecraft  and 
satellite  configurations. 


OPriMAL  ATTITUDE  CONTROL  OF 


AN  ORBITING  SATELLITE 
CONTAINING  FLEXIBLE  APPENDAGES 

I . Introduction 

Background 

In  order  to  meet  functional  requirements,  space  vehicles 
are  frequently  required  to  maintain  a fixed  orientation  with 
respect  to  either  an  orbiting  frame  of  reference  or  some 
inertial  reference  frame.  In  the  early  years  of  space  ex- 
ploration spacecraft  tended  to  be  small,  mechanically  simple, 
and  essentially  rigid.  Today,  however,  satellite  configura- 
tions are  highly  complex  with  many  flexible  devices,  such  as 
antennas , booms , and  solar  panels , which  greatly  affect  the 
attitude  dynamics  of  the  body.  The  stability  and  control 
requirements  for  such  modern  space  systems  are  of  particular 
importance  to  the  designer. 

In  recent  years,  numerous  papers  have  been  presented  to 
determine  the  effect  that  flexible  appendages  have  on  the 
stability  of  space  systems.  Of  particular  interest  are  the 
papers  of  Meirovitch  and  Calico  (Ref  1;  Ref  2),  who  extend 
the  Liapunov  direct  method  of  stability  analysis  to  predict 
the  stability  of  various  spacecraft  configurations.  Their 
basic  strategy  was  to  define  integral  coordinates  to  evaluate 
the  system  Hamiltonian  and  test  for  stability  by  using  the 
Hamiltonian  as  a Liapunov  functional.  Their  research  pro- 
vided a simplified,  alternative  means  of  predicting  stability 
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of  satellites  with  flexible  appendages.  In  the  area  of 
satellite  attitude  control,  methods  of  analysis  are  complex 
requiring  a thorough  investigation  of  a peirticular  vehicle 
configuration  in  order  to  obtain  a prediction  of  optimal 
control  requirements.  By  using  techniques  similar  to  those 
developed  by  Meirovitch  and  Calico  for  predicting  stability, 
this  thesis  will  examine  the  control  requirements  for  satel- 
lites with  flexible  appendages,  and  attempt  to  provide  a new 
method  for  obtaining  the  optimal  control  for  these  space 
systems. 

Problem  Definition 

It  is  well  known  that  a rigid  satellite  is  stable  when 
spun  about  its  axis  of  maximum  moment  of  inertia.  However, 
the  addition  of  flexible  appendages  to  spin  stabilized  space 
systems  can  result  in  an  undesirable  wobbling,  or  coning 
motion,  of  the  vehicle.  To  keep  the  attitude  of  flexible 
satellites  in  proper  orientation  a means  of  controlling  the 
system  is  required.  Frequently,  active  attitude  controllers 
employing  momentum  exchange  or  mass  explus ion  devices  are 
utilized  to  resist  attitude  perturbations.  The  active  con- 
troller can  take  various  forms,  such  as  a momentum  wheel, 
or  as  a reaction  jet  controller  of  the  spin  axis  orienta- 
tion (Ref  3*87).  Regardless  of  what  control  device  is  phys- 
ically implemented,  the  designer  of  flexible  satellites  must 
have  a good  a priori  knowledge  of  the  amount  of  control 
which  will  be  required  to  regulate  a particular  satellite. 

The  overall  direction  of  this  thesis  will  be  to  provide 
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a new,  relatively  simple,  method  for  estimating  the  required 
optimal  control  to  maintain  a particular  spacecraft  attitude. 
The  new  method  will  be  based  on  the  technique  of  integral 
coordinates  previously  employed  by  Meirovitch  and  Calico 
(Ref  1)  to  predict  satellite  stability.  By  using  this  sim- 
plified technique,  the  time  and  effort  required  to  obtain 
an  estimate  for  the  optimal  control  will  be  significantly 
reduced  over  current  complex  methods  employing  modal  analy- 
sis. Having  a simple  means  of  obtaining  the  control  for 
flexible  satellites  would  be  of  great  benefit  to  the  space- 
craft designer,  since  it  provides  a tool  that  is  relatively 
easy  to  employ  upon  which  a control  decision  could  be  made 
during  preliminary  design. 


General  Problem  Formulation 

An  objective  of  this  thesis  is  to  develop  a new  method 
for  obtaining  the  optimal  control  of  a satellite  having 
flexible  appendages.  The  particular  satellite  to  be  inves- 
tigated consists  of  a symmetrical  rigid  body  with  two  flex- 
ible antennas  extending  in  the  positive  and  negative  direc- 
tions of  the  z axis  as  depicted  in  Figure  1.  Mathematical 
analysis  is  performed  on  a discrete  model  (Figure  2a)  and  a 
continuous  model  (Figure  2b)  of  this  satellite  configuration. 
The  discrete  model  will  be  examined  using  spring-mass-damper 
oscillators  to  discretize  the  system  and  the  continuous 
model  will  be  investigated  using  both  a modal  analysis  ap- 
proach and  a simplified  method  of  integral  coordinates. 

The  kinetic  and  potential  energies  will  be  derived  for  the 
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discrete  and  continuous  models  of  this  system.  After  deriv- 
ing these  energy  expressions,  the  system  Hamiltonian  and 
Lagrangian  will  be  formed  and  utilized  in  obtaining  the 
satellite  equations  of  motion.  Once  having  the  equations 
of  motion,  control  will  be  introduced  into  the  problem  along 
with  the  constraint  that  the  satellite  spin  axis  will  main- 
tain a fixed  orientation  with  respect  to  inertial  space. 

To  ascertain  the  amount  of  control  required  to  keep  a nulled 
satellite  position,  an  optimal  control  problem  will  be  for- 
mulated using  discrete,  modal,  and  integral  coordinate 
methods.  The  results  obtained  for  the  three  methods  will 
be  compared  and  an  analysis  will  be  made  regarding  the 
effectiveness  of  using  integral  coordinate  techniques  for 
determing  optimal  control  of  satellites  containing  flexible 
appendages . 

Assumptions 

This  thesis  will  examine  high  spin  stabilized  satellites 
in  earth  orbit.  Since  the  system  is  spinning  at  a relatively 
high  rate,  it  can  be  shown  that  the  torque  due  to  gravity  is 
quite  small.  For  this  reason,  it  will  be  assumed  that  the 
flexible  satellite  under  consideration  is  moment-free  for 
short  time  periods.  This  assumption  is  valid  since  the  ac- 
tual gravity  torque  on  the  system  would  have  an  insignificant 
effect  on  the  attitude  dynamics  of  the  satellite  as  compared 
to  the  effects  associated  with  such  things  as  the  system 
mass  distribution,  system  elastic  properties  and  dynamic 
coupling  thereto,  relative  motion  of  internal  mechanisms. 
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rate  of  conversion  of  mechanical  energy  into  heat  by  dissi- 
pative mechanisms , and  external  forces  of  various  forms 
(Ref  4:1597).  For  this  investigation,  it  will  be  assumed 
that  the  satellite  is  spin  stabilized  about  the  z axis,  and  • 
that  the  mass  distribution,  elastic  properties  of  the  system, 
and  conversion  of  mechanical  energy  into  heat  by  dissipative 
mechanisms  contribute  the  dominant  perturbations  in  the 
attitude  dynamics  of  the  body.  It  will  also  be  assumed  that 
the  angular  momentum  of  the  relative  motion  of  the  satellite 
to  the  earth  is  constant  and  neglegible  compared  to  the 
angular  momentum  of  the  rigid  body  motion.  This  assumption 
is  reasonable  since  this  study  examines  high  spin  stabilized 
satellites.  For  simplicity,  the  center  of  mass  of  the  sat- 
ellite will  be  assumed  to  move  in  a circular  orbit.  Since 
the  gravity  potential  is  approximately  a constant  in  a 
circular  orbit  it  will  not  affect  the  attitude  dynamics  of 
the  body  over  short  time  periods.  Therefore,  the  potential 
energy  will  consist  entirely  of  elastic  strain  energy  of  the 
flexing  antennas.  For  ease  of  calculation,  identical  uniform 
antennas  are  assumed.  Antisymmetric  motion  of  the  satellite 
antennas  will  be  assumed  to  simplify  the  expression  for  the 
kinetic  energy.  It  can  be  shown  that  the  antisymmetric 
motion  assumption  would  represent  the  worst  possible  case 
in  regards  to  the  stability  of  the  flexible  satellite  system 
in  question  (Ref  5‘202).  Additionally,  the  antisymmetric 
motion  assumption  implies  that  the  center  of  mass  remains 
fixed  relative  to  the  main  rigid  body  (Ref  6:1535).  To 
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further  simplify  the  kinetic  energy  expression,  it  will  be 
assumed  that  the  satellite  will  undergo  small  perturbations 
from  the  null  position.  This  is  a valid  assumption  since 
it  can  be  shown  that  for  an  inherently  stable  configuration 
wobbling  of  the  satellite  due  to  perturbations  will  be  con- 
tained in  a small  region  near  the  equilibrium  position.  For 
an  inherently  unstable  configuration  with  no  controller, 
any  perturbation  will  cause  the  cone  angle  of  the  wobble  to 
increase  until  the  body  begins  to  tumble  in  space  (Ref  4 i 1598). 
Since  small  angular  and  vibrational  displacement  perturba- 
tions are  assumed,  angular  and  displacement  terms  in  the 
kinetic  energy  expression  higher  than  second  order  will  be 
ignored.  The  control  required  to  maintain  the  satellite 
system  in  a desired  orientation  will  depend  upon  several 
state  measurements.  It  will  be  assumed  that  any  state  dis- 
placement required  to  implement  the  control  feedback  gains 
can  be  measured  or  somehow  obtained  during  actual  operations 
of  the  satellite.  The  modal  investigation  will  be  assumed 
to  perfectly  model  the  actual  satellite  system,  and  will 
provide  the  reference  for  evaluating  the  effectiveness  of 
the  new  integral  coordinate  technique. 

Sequence  of  Presentation 

The  remainder  of  this  report  is  organized  as  follows. 
Chapter  II  presents  the  mathematical  formulation  of  the 
satellite  system.  The  equations  of  motion  are  derived  for 
the  discrete,  modal,  and  integral  coordinate  techniques  of 
the  uncontrolled  satellite  systems.  The  equations  of  motion 
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are  checked  by  performing  an  eigenvalue  stability  analysis 
and  comparing  the  results  with  previous  investigations  on 
satellites  of  the  same  configuration.  In  Chapter  III, 
control  is  introduced  into  the  system  and  the  optimal  con- 
trol problem  is  formulated.  Closing  this  chapter  is  a brief 
discussion  of  how  control  feedback  gains  may  be  obtained. 

In  Chapter  IV  a presentation  of  the  control  results  is  given. 
The  effectiveness  of  the  integral  coordinate  technique  for 
obtaining  an  estimation  of  satellite  optimal  control  require- 
ments is  analyzed.  The  conclusions  to  be  drawn  from  this 
investigation  are  presented  in  Chapter  V.  A brief  summary 
of  the  investigation  is  presented  in  Chapter  VI.  Also  in- 
cluded in  this  final  chapter  are  recommendations  for  uti- 
lizing integral  coordinate  techniques  for  satellites  with 
flexible  appendages. 


( 

] 

i 


8 


II . Mathematical  Forinulation 


Satellite  in  Earth  Orbit 

As  previously  indicated  the  specific  satellite  config- 
uration under  investigation  consists  of  a symmetrical  rigid 
body  with  two  flexible  appendages  extending  out  the  negative 
and  positive  z axes  (Figure  1).  When  the  satellite  is  in 
the  undeformed  state,  the  x,  y,  and  z body  fixed  axes  rep- 
resent the  principle  moment-of-inertia  axes.  The  mass  mo- 
ments of  inertia  for  the  symmetrical  rigid  body  section  are 
A,  B=A,  and  C about  the  x,  y,  and  z axes  respectively.  The 
satellite  orientation  in  earth  orbit  is  depicted  in  Figure  3* 


i 
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Figure  3*  Satellite  in  Earth  Orbit 
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The  center  of  the  earth  is  considered  an  inertial  reference 
point  with  axes  , and  Z^.  The  radius  vector  RC,  which 

is  constant  in  magnitude  for  a circular  orbit,  is  measured 
from  the  center  of  the  earth  to  the  center  of  mass  of  the 
satellite . 

Discrete  Method  - Equations  of  Motion 

The  discrete  model  of  the  satellite  is  illustrated  in 
Figure  2a.  The  model  consists  of  a symmetrical  rigid  body 
with  spring-mass -damper  oscillators  simulating  the  flexible 
antennas.  The  oscillators  each  have  a mass  m and  are  sym- 
metrically placed  a distance  a along  the  z axis.  The  motion 
of  the  two  oscillators  is  denoted  by  u^,  v^ , U2,  and  V2,  and 
the  angular  velocity  components  of  the  X,Y,Z  system  by 
i'ly,  and  i\.  The  complete  motion  of  the  system  is  described 
by  ten  generalized  coordinates:  three  coordinates  for  the 
motion  of  the  center  of  mass  of  the  satellite,  three  coor- 
dinates for  the  angular  orientation  of  the  body,  and  four 
coordinates  for  the  position  of  the  two  oscillators  with 
respect  to  the  x,y,z  system.  Under  the  assumption  that  the 
center  of  mass  of  the  satellite  moves  around  the  earth  in 
a circular  orbit  and  at  a constant  angular  velocity,  the 
degrees  of  freedom  of  the  system  can  be  reduced  by  three. 

The  assumption  that  the  antennas  move  in  antisymmetric  motion 
implies  that  ui  = -U2  = u and  v^  = -V2  = v.  This  assumption 
allows  a further  reduction  in  the  number  of  generalized 
coordinates  by  two.  Before  the  equations  of  motion  for 
these  generalized  coordinates  can  be  derived,  expressions 
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for  the  kinetic  and  potential  energies  must  be  established 
The  general  expression  of  the  kinetic  energy  for  the 


satellite  system  can  be  written  as 

T = ij  (V*V)dm  (1) 

where  V is  the  inertial  velocity  of  an  element  of  mass  dm. 
For  this  discrete  model,  V can  be  written  in  component  terms 
for  the  rigid  body  and  mass  oscillators  and  the  integral 
evaluated  over  their  respective  domains.  The  general  ex- 
pression for  the  kinetic  energy  becomes 

T = i A'al  + i B'Qy  + i cal  + m(u^  + v^)  + m[v2a| 

+ u^y  + (u^  + v^)az  - 2a^x  2auay  - 2(uv  - 
uv)a2  - Znvilyfly  - 2avaya2  - 2aiaaxPz]  + Const.  (2) 

where  A*  = A + 2ma^,  B'  = B + 2ma^,  and  C are  the  mass  mo- 
ments of  inertia  of  the  undeformed  body  about  the  x,  y,  and 
z axes  respectively.  The  derivation  of  equation  (2)  is 
given  in  Appendix  A. 

Since  the  effect  of  gravity  forces  is  assumed  negli- 
gible, the  potential  energy  consists  entirely  of  elastic 
strain  energy.  For  the  discrete  system  model  the  potential 
energy  expression  is 

V = K(u2  + v2)  (3) 

where  K is  the  spring  stiffness  of  the  oscillators  in  both 
the  X and  y directions.  Here,  the  spring  stiffness  coef- 
ficient K was  chosen  to  be  equal  in  both  directions  since 
uniform  antennas  were  assumed. 


i 


5 


XI 


Energy  dissipation  in  the  system  is  due  totally  to 
damping  in  the  antenna  rods.  In  the  discrete  model,  dis- 
sipation is  modelled  by  placing  damping  pods  on  the  oscil- 
lators. The  equation  expressing  the  effects  of  damping  is 
given  by  Rayleigh's  dissipation  function 

D = d(u^  + ^2)  (4) 

where  d is  the  damping  coefficient  in  the  x and  y directions 
Before  proceeding  to  formulate  the  equations  of  motion, 
it  is  desired  to  express  the  orbitally  referenced  angular 
velocities  liy,  and  in  terms  of  Euler  angles  and  rates 
Figure  4 illustrates  the  3-2-1  rotation  used  in  the  trans- 
formation to  Euler  angles. 


In  Appendix  B it  is  shown  that  the  orbitally  referenced 
angular  velocities  are  related  to  Euler  angles  by 
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(1) 


10  - sin  02 

1 

1 

‘^y 

= 

0 cos  01  sin  01  cos  02 

02 

*^z 

0 - sin  01  cos  01  cos  02 

.*3. 

(5) 


Since  small  angular  displacements  in  and  02  about  a zero 
angle  equilibrium  point  was  assumed,  equation  (5)  can  be 
linearized  by  using  the  Taylor  series  expansion  of  the  sine 
and  cosine.  For  small  angles  Taylor  series  expansion  yields 

sin  01,2  = 01,2 


and 


cos  01,2  = 1 - 1.2 

2 

Substituting  relations  (6)  into  equation  (5)  yields 
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where  0i  and  02  are  angular  velocities  associated  with  small 
perturbation  angles,  and  0^  is  equal  to  the  initial  spin 
rate  about  the  z axis  plus  any  small  perturbation  in  spin 
about  this  axis.  If  equation  (7)  is  substituted  into  the 
kinetic  energy  equation  (2)  and  the  assumption  of  small 
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state  perturbations  is  applied,  terms  of  higher  than  second 
order  in  the  resulting  kinetic  energy  expression  can  be 
ignored,  yielding 

T = i A'9i^  + i B'02^  + i 063^  + m(u^  + v^)  - Bmave^ 

+ 2mau92“  A'il020l  + B'il0i02  + C(i^03  - il0i92) 

+ 2maii92V  + 2ma4]0iu  - 2mk^(uv  - uv)  - 2mai2v02 

- 2nia*2u0i  + i AM^02^  + i B*ii^9l^  + i C(ii^  - 

- 02^^^^)  nul^(u^  + v^)  - 2ma4i^v0i  + 2mai2^u02  (8) 

Since  the  expressions  for  the  kinetic  energy  and  the 
potential  energy  depend  on  generalised  coordinates  and 
velocities  alone,  and  not  explicitly  on  time,  the  system 
is  considered  to  be  non-natural  (Ref  7^77) > For  such  a 
non-natural  system  it  can  be  shown  that  the  Hamiltonian 
assumes  the  form 

I H = T2  - To  + V (9) 

I where  T2  represents  those  kinetic  energy  terms  which  are 

quadratic  in  the  generalized  velocities,  Tq  represents  kin- 
etic energy  terms  which  are  independent  of  the  generalized 
velocities,  and  V is  the  potential  energy  of  the  satellite 
system  (Ref  7s8^t91)*  Therefore,  the  Hamiltonian  becomes 

■ H = i A‘0i^  + i B'02^  + i 003^  + m(u^  + v^)  - 2mav0i 

I + 2mau02  - 2 A'i'1^02^  - i B'iJ^01^  - i C(*i^  - (i^01^ 

- 11^02^)  - nkl^(u^  + v^)  + 2ma4i^v0i  - 2mai‘i^u02 

+ K(u^  + v^)  (10) 


Also  required  for  the  formulation  of  the  equations  of 
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motion  is  the  Lagrangian  which  is  defined  by 

L = T - V (11) 

For  this  particular  system  the  Lagrangian  is  given  by 

L = i A’9i^  + i B'92^  + i C93^  + m(u^  + v^)  - 2raav9i 
+ 2mau92  - A'J929i  + B'ii9l92  + C(ii93  - ii9l92) 

+ 2ma*]92v  + 2mail9iu  - 2mi2(uv  - uv)  - 2mau2v92 

- 2ma.2u9i  + i J B'ii^9i^  + C(ii^  - 91^;^^ 

- 92^^)^)  + mkl  (u^  + v^)  - 2maul^v9i  + 2maii^u92 

- K(u^  + v^)  (12) 

The  system  equations  of  motion  can  now  be  formulated 
by  either  applying  Lagrange's  equations,  which  yield  five 
second  order  differential  equations,  or  Hamilton's  equa- 
tions, which  yield  ten  first  order  differential  equations, 
Hamilton's  equations  are  used  in  this  thesis  to  provide 
five  first  order  differential  equations  of  the  generalized 
coordinates  (9i,  92,  63,  u,  and  v)  and  five  first  order 
differential  equations  of  the  generalized  momentum  (P9i, 

P92*  P931  Pu*  3.nd  Pv)-  The  generalized  momentum  is  defined 

Pi  s iL_  i = l,"-.5  (13) 

s^i 

For  a holonomic  system  subject  to  forces  not  derivable  from 
a potential  function,  Hamilton's  equations  of  motion  take 
on  a special  form.  In  particular,  when  the  nonconservative 
forces  in  the  system  are  derivable  from  Rayleigh's  dissipa- 
tion function,  Hamilton's  equations  become 
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i = 1.  — .5  (14) 


qi  = 3H_ 


- ^ 

3qi  3qi 


where  Rayleigh's  dissipation  function  D is  expressed  in 
equation  (4)  (Ref  7«95)*  A detailed  derivation  of  the  ten 
first  order  differential  equations  is  presented  in  Appendix  C. 
These  equations  of  motion  are  summari:ied  below i 


01  = Pei  + aPv  + iiep 

A 

02  = P92  - aPu  + (c  - l)il0i 

B B 


e 


3 = -Q 

C 


^ _ aP02  + B'Pu  + ilv  - Caflei 

B 2mB  B 


V = fPn  + A^v  - Qu 
A 2mA 


P0  1 = (1  - C)ca  01  + (1  - C)ilP02  + CailPu 

^ B B B 


P02  ~ ” Cil  02  “ HP 01 


P93  = 0 


Pu  = - 2Ku  - 2Du 


Pv  = - ilPu  - 2Kv  - 2Dv 


(15) 


Note  that  equation  (15)  can  be  expressed  in  state  notation  as 


x(t)  = Ax(t)  (16) 

where  x(t)  is  a 10  x 1 state  vector 

A is  a constant  coefficient  10  x 10  matrix 
A stability  analysis  of  the  discrete  system  can  now  be 
accomplished  by  evaluating  the  eigenvalues  and  eigenvectors 
of  the  A matrix.  This  task  was  performed  in  this  thesis  to 
check  the  validity  of  the  equations  of  motion  shown  in  equa- 
tion (15)-  Validity  was  confirmed  by  comparing  stability 
results  numerically  calculated  in  this  thesis  to  the  sta- 
bility results  obtained  by  Meirovitch  in  Ref.  4 for  the  same 
satellite  system.  Although  Ref.  4 used  Lagrange's  equations 
to  formulate  the  equations  of  motion,  the  stability  results 
for  various  system  parameters  were  identical  to  the  stabil- 
ity results  obtained  by  using  Hamilton's  equations.  Fig- 
ures (5)  and  (6)  depict  ratio  of  spin  rates  obtained  from 
Ref.  4 and  from  this  thesis  for  particular  satellite  para- 
meter sets.  Additionally,  eigenvalue  stability  analysis  was 
in  close  agreement  with  that  obtained  in  Ref  4.  Figure  7 
illustrates  the  stability  regions  obtained  in  Ref.  4 and 
in  this  thesis.  Note  that  for  these  stability  plots,  stable 
regions  are  under  the  curve  while  unstable  regions  are  above 
the  curve.  iJhile  not  meant  to  be  an  exact  proof,  the  simi- 
larity of  stability  results  obtained  in  this  analysis  tends 
to  support  the  validity  of  the  equations  of  motion  given 
by  equation  ( 15) • 
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Figure  5.  Stable  System  Response:  a)  (Ref  4);  b)  Thesis 
Investigation 
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Figure  ?.  Stability  Plot  - Discrete  System:  a)  (Hef  4); 
b)  Thesis  Investigation 
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From  this  point  in  the  investigation,  some  type  of 
control  could  be  introduced  into  the  discrete  system  to 


keep  a perturbed  satellite  in  a proper  or  desired  orienta- 
tion. However,  before  proceeding  with  the  control  inves- 
tigation, the  equations  of  motion  of  both  the  modal  analysis 
and  integral  coordinate  techniques  will  be  derived.  After 
obtaining  these  attitude  motion  equations,  control  will  be 
introduced  and  the  results  of  each  method  will  be  presented 
and  compared. 

Continuous  Modal  Method  - Equations  of  Motion 

The  continuous  model  of  the  satellite  system  is  illus- 
trated in  Figure  2b.  This  model  consists  of  two  flexible 
uniform  rods  of  length  1 connected  to  a symmetric  rigid 
main  body.  The  rods  are  attached  at  a distance  h from  the 
center  of  mass  of  the  system.  Each  antenna  rod  has  a mass 
per  unit  length  of  p and  a total  mass  of  m = p x 1.  Pro- 
ceeding as  in  the  discrete  system  analysis,  the  kinetic  and 
potential  energy  expressions  will  now  be  formulated. 

The  general  expression  for  the  kinetic  energy  is  given 
by  equation  (1).  Using  the  same  assumptions  as  in  the  dis- 
crete system  analysis,  the  kinetic  energy  for  the  continuous 
system  can  be  written  as 

T = i + t [u^  + v^ 

+ + (u^  + v^).ii,2  - 2zviij^  + Pziiiiy 

- 2(uv  - uv)*ij.  - 2uVi2x»iy  - 2zVi2yii^  - 2zu*2x*lz]  pdz 

+ Const.  (17) 
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In  equation  (1?)  A*  and  B'  are  the  mass  moments  of  inertia 
about  the  x and  y axes  and  includes  the  mass  of  the  rods  in 
the  undeformed  position.  If  the  transformation  presented  on 
page  13  is  again  performed,  the  kinetic  energy  can  be  writ- 
ten in  terms  of  Euler  angles  and  rates.  If  the  assumption 
of  small  state  perturbations  is  applied,  terms  of  higher 
than  second  order  can  be  ignored,  resulting  in  the  following 
expression  for  the  kinetic  energy 

T = i A'(0i^  - 29291^1  + 92^1^)  + i B'(e2^  + 

+ + i 0(93^  - 29i92ii  + 293,1  - 9i^a^  - 

920^  + ii^)  + I ^ + (u^  + - 

^ h 

2zv(9i  - 92ii)  + 2zu(92  + 9iii)  - 2(uv  - uv)!!  - 
2zv(92ii  + 0ii]^)  - 2zu(9j^Q  - 92ii^)J  pdz  + Const.  (18) 


Breaking  the  kinetic  energy  into  T2  terms  (quadratic  in 
generalized  velocities),  terms  (linear  in  generalized 
velocities),  and  Tq  terms  (containing  generalized  coordin- 


ates alone)  yields 


h + 1 


T2  = i A'9i^  + i B'9'2^  + 2 C9*3^  + \ [u^  + v^ 

h 

- 2zv9j^  + 2zu92]pdz 

Pi  = - A'.i929i  + B',i9i02  + C(Ai93  - ,19102) 


,h  + 1 


[2z,192V  + 2Zii9iu  - 2,l(uv  - uv) 

- ~ 2Zilu9iJpd. 
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'o  = i ^ + i - 01^  - 02^) 


h + 1, 


+ j ■*■  “ 2z*2^v0j^  + 2zil^u023pdz  (19) 

For  the  continuous  system  the  potential  energy  expres- 
sion, which  is  considered  to  be  caused  only  by  the  elasti- 


city of  the  rods , is 


V = f ^ ^£I  /9£u\^  + /i2v\^  dz  (20: 

^h  L!^8z2;  \^9z2;_ 

where  it  is  assumed  that  the  bending  stiffness  of  the  rods 
are  identical  in  both  the  x and  y directions  (Ref  8j229). 

The  damping  in  the  system  is  again  given  by  Rayleigh's 
dissipation  function.  For  the  continuous  system  the  dis- 
sipation function  has  the  form 

D = / d(u^  + v'^)dz  (21] 


where  d is  the  damping  coefficient  in  the  x and  y directions. 

The  general  form  of  the  Hamiltonian  and  the  Lagrangian 
remain  as  given  by  equations  (9)  and  (11).  For  the  contin- 
uous system  the  Hamiltonian  can  be  written  as 

H = i A'0i^  + i B'02^  + i C03^ 

- 2zv0i  + 2zu02]pdz  - i A'il^02^  “ i B*ii^01^ 

+ i CL2^(ei^  + 02^  - 1)  - ^ ^[(u^  + v^)a^ 

h 

- 2zil^v0i  + 2zi)^u02lpdz  + f El  f 3 ^u\ 

■'  h \dz2/ 

+ f iZl]  dz  (22) 

[3z2/_ 
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The  corresponding  Lagrangian  can  be  written  as 

L = i A'9i^  + i B'92^  + i C93^  + f [u 

■'  h 


2 ^ .2 
+ V 


- Zzvk'i  + 2zu923pd*.  - A'ii929i  + B*ii9l92  + 
.h  + 1, 


Cii(93  - 9i92)  + j L25iii92V  + 2zii9iu  - 


2i2{uv  - uv)  - 2ziiv92  - 2ziiu9jpdz  + ^ A’ii^92^ 


h + 1, 


+ i B'U^9i^  + i Ca^d  - 9i^  - 92^)  + f [(u 

h 

+ v^)i2^  - 2zii2^v9i  + 2zii^u923pdz  - 


h + 1 

I El 

^ h 

V9^2;  1 

[^z2)  J 

dz  (23) 


It  can  be  noted  that  in  equations  (22)  and  (23)  elastic 
displacements  appear  in  integrals  defined  over  the  elastic 
domain.  The  spatial  and  time  dependency  of  u and  v presents 
some  difficulty  in  analyzing  the  continuous  system  in  that 
the  Hamiltonian  and  Lagrangian  contain  both  continuous  and 
discrete  coordinates  as  well  as  spatial  derivatives  of  the 
continuous  coordinates.  Such  systems  are  termed  hybrid, 
that  is,  they  give  rise  to  both  ordinary  and  partial  differ- 
ential equations  in  describing  the  motion.  In  order  to  cir- 
cumvent this,  the  system  will  be  discretized  by  approximat- 
ing elastic  displacements  by  means  of  modal  analysis. 

In  this  thesis  the  first  four  natural  modes  of  antenna 
motion  are  investigated.  Using  the  antisymmetric  motion  as- 
sumption, the  modal  amplitudes  in  the  x and  y directions  can 
be  described  by  u^,  U2,  U3,  U4,  vj^ , V2,  V3,  and  V4.  Using 
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the  assumption  that  attitude  motion  need  only  be  considered, 
the  complete  attitude  orientation  of  the  satellite  can  be 
described  by  eight  natural  mode  displacements  and  three 
angular  rotations  of  the  body.  The  modal  analysis  treats 
the  elastic  antennas  as  flexible  beams  subject  to  the  bound- 
ary condition  that  the  beams  are  rigidly  attached  to  the 
satellite  at  one  end  and  free  at  the  other.  The  elastic 
motion  of  the  individual  antenna  can  be  considered  to  con- 
sist of  a superposition  of  the  normal  modes  of  a beam  having 
the  same  mass  and  stiffness  distribution  as  the  flexible 
antenna  (Ref  4:l601).  In  terms  of  an  equation,  the  overall 
motion  in  the  u and  v directions  is  described  by 

n 

U(2,t)  = z $i(2)ui(t) 

i = 1 

n 

V(2,t)  = Z $i(2)vi(t)  (24) 

i = 1 

where  z represents  any  point  along  the  beam.  Here,  the 
normal  modes  or  eigenfunctions  associated  with  a uniform 
cantilever  beam  satisfy  the  differential  equation 

El  ^^^i  - = 0 i = l, ,n  (25) 

where  the  effects  of  shear  deformations  and  rotary  inertia 
are  neglected  (Ref  9!l98)«  For  a uniform  beam  the  quanti- 
ties £I  and  p are  constants.  In  equation  (25)  (jj^  is  the 
natural  frequency  of  the  i"^^  mode.  Additionally,  the 
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boundary  conditions  on  equation  (25)  for  a clamped-free  beam 
translate  to 

= 0,  ^^i  = 0 z = h 

dz 

2 2 

El  ^ = 0,  d(EI  ^ ^i)  /dz  = 0 z = h+  l (26) 

dz2  dz^ 

previous  solutions  to  the  problem  associated  with  equations 
(25)  and  (26)  yield  eigenfunctions 

<E>i  = [sin  3il  - sinh  Bjl]  [sin  Bi(z  - h)  - sinh  ei(z  - h)] 

V pi  sin  6il  sinh  6il 

[cos  Bj^l  + cosh  Bj^l]  [cos  ^^(z  - h)  - cosh  3i(z  - 

V pi  sin  6il  sinh  3 il 

in  which 

(BiD^  = (Di^pl  , (28) 

El 

where  3^1  is  determined  by 

cos  3il  cosh  3il  = -1  (29) 

Values  of  3j^l  for  various  beams  are  readily  obtainable  from 
reference  text  books  in  structures.  For  the  clamped-free 
beam  these  values  are 

3il  = 1.875 

62I  = 4.694 
63I  = 7.855 

3il  = 2i  - 1-  (30) 

2 “ 
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In  addition,  the  eigenfunctions  of  equation  (2?)  are  ortho- 
gonal and  can  be  normalized  with  respect  to  the  mass  per 
unit  length  function  p such  that 

1 

J p<fi(  z)<|>j(  z)dz  = 5ij  (31, 

where  represents  the  Kronecker  delta  (Ref  li35)« 

Using  the  above  modal  analysis  techniques,  integral 
terms  in  the  Hamiltonian  and  Lagrangian  can  be  discretized. 


For  example, 


.h  + 1 


f u^pdz  = s Ui' 
-'h  i = 1 


^-2  . ” • 

I u'^pdz  = Z Ui 

Jk  i = 1 


h + 1. 

I livpdz 

■’  h 

h + 1 

I upzdz 

^ h 

j 2z92updz 

J h 


Z UiVi 
i = 1 


2 ^ 

(ml'^)  z Sz^ui 
i = 1 

5 i n 

202(ml'^)“  Z Sz^ui 
i = 1 


where 


^ (Sil)^sin  Bfl  sinh  Bil  ~ 

(Ref  li5^)  (33) 


Additionally,  it  can  be  shown  (Ref  8j220)  that 


p 2 ^ o o ^ o 

El  + /l£v\  dz  = Z ooi^ui^  + Z 

h LUz2j  (sz2j  J i = 1 i = 1 


i^  (34) 
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Since  higher  order  modes  will  contribute  minimally  to 
the  overall  antenna  motion,  this  investigation  will  restrict 
analysis  to  the  first  four  natural  modes  of  motion.  With 


this  in  mind,  the  expression  for  the  Hamiltonian  becomes 

4 4 

H = i A’9i2  + I B'02^  + i 003^  + Z ui^  + z 

i = 1 i = 1 

- 20i(ml^)  Z Szivi  + 202(ml^)  Z SziUi 

i = 1 i = 1 

- i + i c,22(9^2  + q^2  _ 

- ^ ui^  - Z vi^  + 20lil^(ml^)^  Z SziVi 

i=l  i=l  i = l 

- Z SziUi  + Z (bi^(ui^  + vi^)  (35) 

i = 1 i = 1 


Similarly  the  Lagrangian  becomes 

L = i A'0i^  + i B’02^  + i 003^  + Z iui^  + vi^) 

i = 1 

1 4 

+ 2(ml^)^  Z Szi(ui02  - Vi0i)  - A’a02ei  + BM0i02 
i = 1 


4 

+ 011(03  ~ ®1®2)  2i2(ml^)®  Z S2i(v202  + Ui0i  - 

i = 1 
4 

Vi92  - Ui0i  ) + 2il  Z (uivi  - uivi)  + ^ A'il^02^  + 

i = 1 

i B’il^9l^  + i Cil^d  - 01^  - 02^)  -t-  ^ (ui^  + Vi^) 

i = 1 

i 

+ 2i2^(ml^)^  i;  S2.(ui02  - Vi0i)  - cui^  z (ui^  + Vi^) 
i = 1 i = 1 

(36) 


One  is  now  in  a position  to  obtain  the  equations  of 
motion  for  the  continuous  system  by  the  modal  method.  The 


generalized  momentum  is  again  given  by 

Pk  = -SL  ^ = 1,  — ,11  (37) 

^‘ik 


With  Rayleigh's  dissipation  function  D given  by  equation 
(21),  Hamilton's  equations  become 

4k  = _M  • Pk  = - ^ = 1,  — ,11  (38) 

3Pk  S'lk  S^lk 


The  resulting  22  first  order  differential  equations  of 
motion  are  derived  in  Appendix  D and  summarized  below 

01  = 02  + PiL 

Da 

1 4 

02  = (_C  - 1)  9i  f £02  - (ml£)^  Z SziPui 

Db  Db  Db  i = 1 

Vi  = - ilu-i  + (ml^)  S2iP0.  + Szi  ^ SziPvi 

J ''  Da  ^ ^ ^ ^i  = 1 

+ j = 1.  — .4 

2 

ui  = - (!!i1£)^Sz.P02  - (ml£)^CaSz.9i  + 

•'  Db  Db  ^ 

2 ^ 

+ s I s p j = 1,  — ,4 

9 2 i ^ 

Pei  = (1  - _c)c*rei  + (i  - _c)i]P92  + (nii_)  Cii  z Sz.Pui 

Db  Db  Db  i = 1 

P92  = " Ci1^92  - 
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y 


2u)i2v:i  + gdilu;.  - 2d  (ml£)^S2:iP9i-  ilPui 

J p p Da  J ^ j 

|[Pvj  + 2ra^  Sz-  ^ I ^Sz^PvJ  j = 1.— 


.4 


^uo 


- 2(i)T^UT  - ^ ^ (inl_)  3z-(Ci^9l  + Pq-)  + ilPyi 

p p Db  •’ 


■ d r Pu-;  + 2ml£  3z^  I SziPu^l  j = 1, 

p L J Db  •'  i = 1 “■ 


.4 


P 

>3  = Ih.  - 


Pe^  = 0 (39) 

The  relations  appearing  in  equation  (39)  were  verified 
in  much  the  same  manner  as  the  discrete  system  equations. 
Noting  that  O3  and  P93  uncouple  from  the  system,  an  eigen- 
value stability  analysis  was  performed  on  the  remaining 
20  X 20  constant  coefficient  matrix.  Again,  the  results 
obtained  in  this  analysis  correspond  to  the  modal  method 
results  obtained  in  Ref.  4 for  the  continuous  system.  This 
correspondence  of  results  is  illustrated  in  Figure  8 where 
similar  stability  plots  are  depicted  for  various  system 
parameters.  In  Figure  8,  stable  regions  are  below  the  curve 
and  unstable  regions  are  above  the  curve. 

For  this  investigation,  the  equations  of  motion  ob- 
tained by  the  above  modal  technique  will  serve  as  a true 
model  of  the  satellite  system.  Later  in  this  thesis,-  con- 
trol will  be  added  to  the  system  and  results  obtained  using 
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Figure  8.  Stability  Plot  - Continuous  System  a)  (Ref  4) 
b)  Thesis  Investigation 
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integral  coordinate  equations  of  motion,  which  are  derived 
in  the  next  section,  are  compared  to  results  obtained  by 
this  so-called  true  model. 


Integral  Coordinate  Method  - Equations  of  Motion 

As  can  be  noted,  the  modal  analysis  method  of  the  pre- 
vious section  leads  to  a rather  large  and  complex  system 
of  differential  equations.  Additionally,  modal  analysis 
does  leave  some  question  as  to  the  accuracy  of  results 
since  series  truncation  is  utilized.  In  this  section,  a 
method  of  integral  coordinates  will  be  introduced  which  will 
significantly  reduce  the  number  of  differential  equations 
to  be  solved.  Although  the  integral  coordinate  technique 
has  fewer  differential  equations,  the  system  is  still  based 
on  a continuous  representation  with  results  maintaining  a 
high  degree  of  correctness. 

The  model  for  the  satellite  system  remains  that  given 
by  Figure  2b.  Since  the  model  is  continuous,  the  general 
expressions  for  the  system  Hamiltonian  and  Lagrangian  are 
once  again  given  by  equations  (22)  and  (23).  In  order  to 
circumvent  the  problems  associated  with  evaluating  the  spa- 
tially dependent  integrals  in  equations  (22)  and  (23),  new 
coordinates  will  be  selected.  To  this  end,  the  following 
integral  coordinates  are  defined 


h + 1 

u = I puzdz  , 

^ h 

h + 1 

V = I pvzdz  , 

h 


u = / puzdz 

•'  h 

V = I pvzdz  (^0) 

^h 
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For  the  integral  coordinate  formulation,  it  can  be  surmised 
that  the  general  coordinates  will  be  02,  63,  u,  and  v. 
The  immediate  advantage  of  this  approach  is  in  its  simplic- 
ity. The  continuous  system  may  now  be  investigated  using 
the  same  number  of  generalized  coordinates  as  in  the  dis- 
crete system,  while  maintaining  the  completeness  of  the  more 
complex  modal  analysis  investigation. 

In  conjunction  with  the  definitions  given  in  equation 
(40),  integral  terms  in  equations  (22)  and  (23)  which  are 
quadratic  in  u or  v may  be  handled  by  using  Schwarz's  in- 
equality. This  inequality  relationship  is 

M pvzdzj  £ J pz'^dz  J pv^dz  (Ref  1*59)  (41) 
or 


^ o / ^ ^ \2 

I pv^dz  5 ( r pvzdzj 

a.  1 


^ 2 


(42) 


Approximating  by  using  the  equality  sign  above,  and  recall- 
ing the  definition  of  v yields 

.h  + 1 


L 


pv  dz  2! 


_ 2 
V 


(43) 


/ 


h + 1 


2 

pz  dz 


Similarly, 


/ 


h + 1 


pu^dz 


_ 2 
u 


/ 


h + 1 


pz^dz 
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(44) 


W 


If  the  same  procedure  is  followed  for  integrals  with  qua- 
dratic terms  in  u and  v,  substitution  into  equation  (22) 
for  the  Hamiltonian  yields 


-I.*?  ^ *2  ^ * 2 *2  *2  ••  •• 

H = i J 3*92  + i CQy  + I_  + vf  - 20iv  + 292U 

RI  RI 


+ (ti)i 


where 


- i B'0i’a^  - i Cil^(l  - 01^  - 02^) 

- ClH 

^ + v£]  + 29i;]^v  - 292a^u 

RI  RiJ 

(45) 

RI  = 

rh  + 1 2 

I pz  dz 

J h 

A'  = 

A + 2m(h  + 1)^  + 2ml2 

2 12 

B*  = 

A* 

(46) 

Similarly,  substitution  into  equation  (23)  for  the 
Lagrangian  yields  the  following 


L = i A’ei^  + i B'92^  + i C03^  + if  + if  - 201^  + 202U 

RI  RI 


- A *02011^  + B*9i92i')  - 00102^'^  + 00311  + 202*lv  + 29iilu 

- 202^1v  - 20ii'E  + i A'02^il^  + i i 00^(1  - 


01^  - 02^)  + (ii^  “ uf  + vfl  + 202il^u  - 29ii’)^v 

.RI  RiJ 


(^7) 


For  simplicity,  the  term  (uv  - uv)  in  the  Lagrangian  has 
been  intentionally  neglected  due  to  its  small  overall  effect 
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on  the  system.  The  differential  equations  of  motion  for 
the  above  integral  coordinate  technique  can  now  be  formu- 
lated by  again  applying  Hamilton's  equations  given  in  equa- 
tions (13)  and  (14).  With  Rayleigh’s  dissipation  function 
defined  by  equation  (21),  it  is  shown  in  Appendix  E that 
the  integral  coordinate  equations  of  motion  become 


01  = Pei  + RI  Py  ^ -202  + 2ilu 

A*  - 2RI  A*  - 2RI 

02  = ^02  ~ Pu  + ( C - l)il0i  + 2Qv 

B*  - 2RI  B*  - 2RI  B’  - 2RI 

^ = B'RI  P^;-  _ RI  P02  - CilRI  01  _ 2ilRI  V 

2(B’  - 2RI)  B*  - 2RI  B'  - 2RI  B'  - 2RI 

^ ^ A'RI  ?;7  RI  P9i  ^ 2nRI  u 

2(A'  - 2RI)  A*  - 2RI  A'  - 2RI 


(1  - c )ct')  ©1  + (1  - c )apflo 

B*  - 2RI  B'  - 2RI 


+ COR  I p^  _ 2ca^7 

B'  - 2RI  B*  - 2RI 


P92  = - OPei  - ca^92 


P~  = 


2a1P01  - 2^)RI 

A'  - 2kI  A'  - 2RI 


4i^^u  + 2(a^ 

A*  - 2RI 


2x  - 

'*>1  ) _li 
RI 


P7  = 


+ d (2Ci)0l  + 2P02  - B'P-  + 4tlv) 

p(B'  - 2RI) 

2nP92  + 2.1^1  Pu  - 2C.)20i  _ ^ 

B'  - 2RI  B*  - 2RI  B'  - 2RI  B*  - 2RI 

2(;1^  - u)i^)  J.  - d (A'P7  - 2P0,  - 4il^r) 

RI  p(B’  - 2RI)  ^ 


03 


^^3  - il 
C 


(48) 


Pe3  = 0 

The  equations  of  motion  are  now  derived  for  the  three 
methods  under  investigation.  In  the  next  chapter,  control 
will  be  introduced  into  the  problem.  After  the  introduction  j 

of  control,  the  control  requirements  for  the  three  methods  j 

i 

will  be  determined  and  the  effectiveness  of  the  integral  | 

coordinate  technique  for  obtaining  control  feedback  gains  j 

will  be  analyzed.  | 
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III . Control 

Up  to  now  the  only  control  on  the  satellite  has  been 
the  passive  attitude  control  inherent  in  spin  stabilized 
systems.  Active  attitude  control  of  the  satellite  is  re- 
quired for  two  reasons.  For  one,  even  though  a stable 
spinning  satellite  will  tend  to  resist  attitude  perturba- 
tions, in  time,  its  spin  axis  will  move  away  from  a desired 
direction  due  to  the  influence  of  a gravity  gradient,  aero- 
dynamic force,  or  other  such  perturbative  force.  Secondly, 
active  attitude  control  is  required  to  reorient  the  spin 
axis  of  the  satellite  if  a change  in  pointing  direction  is 
required  during  flight.  In  this  thesis,  the  satellite  will 
be  initially  displaced  slightly  from  a null  position  to 
simulate  the  effect  of  either  a perturbative  force  or  a 
desired  change  in  spin  axis  orientation.  To  effect  the 
attitude  change,  control  will  be  applied  to  the  angular 
velocities  and  02  of  the  satellite  system.  A momentum 
exchange  or  a mass  explusion  device  could  be  utilized  to 
impart  this  change  in  angular  velocities,  although,  the  ac- 
tual physical  mechanism  is  immaterial  to  this  investigation. 

Optimal  Control 

Modern  control  theory  techniques  will  be  applied  to 
obtain  an  estimation  of  the  optimal  control  required  to  ; 

produce  desired  attitude  changes  in  the  satellite  system.  : 

In  this  thesis,  control  is  considered  optimal  in  the  sense 
that  minimum  angular  velocity  control  is  used  to  minimize 
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state  deviations.  For  all  three  methods  under  investiga- 
tion, the  uncontrolled  satellite  is  described  by  equation 
(l6).  Since  the  system  is  linear  and  time-invariant,  the 
differential  equation  for  the  controlled  satellite  system 
is  given  by  the  matrix  equation 
• 

x(t)  = Ax(t)  + Bu(t)  (Ref  10j148)  (49) 

where 

A is  a (n  X n)  plant  matrix 
B is  a (n  X r)  input  matrix 
x(t)  is  a (n  X 1)  state  vector 
u(t)  is  a (r  X 1)  control  vector 
Optimal  control  is  now  sought  such  that  a quadratic  cost 
functional  involving  the  weighted  components  of  the  states 
plus  the  control  is  minimized  when  the  final  time  is  unspec- 
ified. In  terms  of  an  equation,  the  quadratic  cost  function 
is  given  by 

J = if  (x'^(t)Qx(t)  + u^(  t)Ru(t)  )dt  (50) 

■'o 

where  Q and  R are  positive  definite  symmetric  weighting 
matrices  on  the  states  and  control  vectors  respectively 
(Ref  10:149).  If  4 is  large  it  would  imply  that  the  corre- 
sponding state  components  will  be  rapidly  regulated.  If 
there  are  large  entries  in  R,  this  would  imply  that  energy 
in  the  corresponding  control  element  will  be  kept  small. 

In  this  thesis,  both  Q and  R are  set  equal  to  the  identity 
matrix,  (I.),  implying  equal  weighting  for  the  states  and 
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the  control.  For  the  cost  function  and  state  relationship 
presented  above,  it  can  be  shown  that  a statement  of  the 
optimal  control  is  given  by 

u(t)  = - R"Vpx(t)  (51) 

where  the  determination  of  the  optimal  control  involves  the 
solution  of  P by  the  matrix  Riccati  equation  (Ref  10j151)- 
The  matrix  Riccati  equation  is  given  by  the  following  dif- 
ferential equation 

P + PA  + a'^P  - PBr’Vp  = 0 (52) 

where  the  steady  state  solution  can  be  determined  by  setting 
P equal  to  zero.  Defining  the  feedback  gain  matrix  as 

-1  T 

F = - r bp  (53) 

the  closed-loop  system  matrix  differential  equations  can 
be  written 

x(t)  = (A  + BF)x(t)  (5^) 

For  the  closed-loop  system,  stability  can  be  determined  by 
numerically  calculating  the  eigenvalues  and  eigenvectors 
associated  with  the  time-invariant  matrix  coefficient  on 
x(t)  in  equation  (5^)*  A block  diagram  of  the  closed-loop 
system  is  illustrated  in  Figure  9. 


i 
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Figure  9.  Closed-Loop  System  of  Control  Feedback 

Control  Feedback  Gains 

In  order  to  obtain  control  feedback  gains  for  the 
three  methods  investigated,  values  of  R and  B were  selected 
and  steady  state  solutions  to  the  matrix  Riccati  equation 
were  numerically  determined.  The  AFIT  Subprogram  OPTCON 
(Ref  11)  and  the  Aerospace  Medical  Research  Laboratory 
Subroutine  MRIC  (Ref  12)  were  used  in  obtaining  the  numer- 
ical solutions  to  the  steady  state  Riccati  equation.  With 
the  ability  to  obtain  feedback  gain  matricies  for  particular 
initial  conditions  and  satellite  parameters,  it  is  now 


possible  to  obtain  and  compare  control  requirements  for  the 
three  techniques  being  investigated. 


IV . Results  and  Analysis 


In  chapter  II  the  dynamics  of  the  satellite  were  formu- 
lated for  the  three  techniques  being  investigated.  A means 
of  controlling  the  system  was  introduced  in  chapter  III. 

In  this  chapter,  specific  examples  are  examined  and  the 
results  for  the  various  techniques  are  presented. 

For  each  method  being  investigated,  optimal  control  was 
determined  by  numerically  obtaining  control  feedback  gains, 
F,  based  on  equations  (52)  and  (53)*  As  previously  men- 
tioned, Subroutine  MRIC  and  Subprogram  OPTCON  were  utilized 
to  obtain  feedback  gains  for  the  system.  These  feedback 
gains  were  substituted  into  the  general  closed-loop  system 
differential  equation  (5^)#  which  describes  the  satellite's 
controlled  motion.  The  closed-loop  differential  equation 
was  then  integrated  forward  in  time  by  using  the  CC6600 
Subroutine  ODE  (Ref  13)*  For  the  example  which  follows,  the 
satellite’s  angular  displacements  and  antenna  motion  are 
plotted  versus  time  for  both  uncontrolled  and  controlled 
situations . 

Example  1 

In  this  example  a stable  satellite  configuration  was 
examined.  The  following  initial  conditions,  satellite  pa- 
rameters, and  weighting  matrices  apply  to  all  three  methods. 

Initial  Conditions: 

Ql  = 0° 

02  = 3°  = .0523598776  Radians 
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u 


= 0 ft  = u = ui  = U2  = U3 
V - 0 ft  = V = = V2  = V3 
t = 0 sec 

Satellite  Parameters: 
iig  = 50  Rad/sec 


^n 

= = 60  Rad/s 

II 

00. 

C 

= 60  slugs-ft^ 

C/A 

= 1.5 

B 

= A 

RA 

= .1 

a 

= 2 ft 

h = 2 ft 
1 = 1 ft 

Zeta  = d/2m'a)j^  = .5 
d/p  = 10 

Weighting  Matrices: 

Q = [I] 

R = [I] 


U4 

V4 


Discrete  System  Results 

For  the  discrete  system  model,  Figures  10  and  11 
illustrate  the  satellite's  motion  for  the  uncontrolled 
system.  It  can  be  noted  that  for  an  inherently  stable  sys- 
tem, the  antenna  displacements  of  u and  v tend  to  dampen 
out  over  a period  of  time.  The  damped  motion  of  the  an- 
tenna was  characteristic  for  any  value  of  zeta  greater  than 
zero.  Recalling  that  and  02  are  angular  measurements  of 
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Figure  11.  Uncontrolled  Angular  Motion  - Discrete  System 
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the  body  fixed  axis  with  respect  to  an  inertial  axis,  it 
can  be  seen  from  Figure  11  that  the  spin  axis  assumes  a 
new  pointing  direction.  A controller  is  required  to  drive 
these  angular  displacements  back  to  zero  so  that  the  spin 
axis  is  once  again  oriented  to  the  null  position.  Table  I 
provides  the  feedback  gains  required  to  control  the  system. 


Table  I. 

Discrete  Feedback  Matrix 


Col  1 

Col  2 1 

Col  3 1 

Col  4 j 

Rov.- 

1 

-3.7003r:+02 

-4.5512F+00  -2.7374i-05 

-I.032OS+OI 

Row 

2 

-2.737^2-05 

7.299OE-06  -5.4318E+02 

7.09582-06 

0015 

1 

col  7 

Col  8 

Row 

1 

9.90532-01 

2.26242:-02  -1.2334";-01 

4.75342-01 

Row 

2 

- I.8I06S-OI 

-7.03792-08 

9.83472-01 

-4.45722-07 

The  optimal  control  at  any  instant  of  time  is  just  the  matrix 
multiplication  of  the  feedback  gain,  F,  times  the  current 
state  vector. 

u(t)  = Fx(t)  (55) 

Recalling  that  0-3  and  uncouple  for  all  three  methods 
investigated,  the  states  of  the  discrete  system  can  be  des- 
cribed by  an  (8  x 1)  vector  written  in  the  following  order. 


xi  = 01 
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(56) 


X5  = P01 

X6  = Pv 
X7  = P92 

X8  = Pu 


Figure  12  illustrates  the  angular  displacements  and  antenna 
motion  of  the  controlled  discrete  system.  As  can  be  ob- 
served, the  satellite  is  returned  to  the  desired  null  ori- 
entation in  just  over  one  second.  For  this  particular  set 
of  satellite  parameters,  the  stability  of  the  closed-loop 
system  is  shown  by  the  all  negative  eigenvalues  illustrated 
in  Table  II. 


Table  II. 

Discrete  Closed-Loop  Eigenvalues 


Real 

Imaginary 

-.33957E+01 

-.52787E+OI 

-.33957E+01 

.52787E+OI 

-.5417IE+02 

-.I0598E+03 

-.5417IE+02 

.IO598E+03 

-.19348E+03 

-.I9573E+03 

-.19348E+03 

.I9573E+03 

-.27159E+03 

-.27615E+03 

-.27159E+03 

.276I5E+O3 
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Figure  12.  Coutrolled  States  - Discrete  System 


Modal  Analysis  Results 

The  modal  analysis  investigation  of  the  continuous 
system  model  was  the  second  method  applied  to  this  exaunple 
problem.  In  order  to  obtain  an  insight  into  the  unforced 
satellite  motion,  the  equations  of  motion  were  numerically 
integrated  and  the  uncontrolled  states  were  plotted  against 
time.  Figure  13  illustrates  the  uncontrolled  angular  dis- 
placements of  0]_  and  0 2.  A breakdown  of  the  uncontrolled 
motion  for  the  first  four  modes  of  antenna  vibration  is 
presented  in  Figures  14  through  1?.  As  expected,  the  first 
mode  significantly  dominates  the  antenna  motion.  The  magni- 
tude of  the  overall  uncontrolled  modal  displacements  in  the 
u and  V directions  is  illustrated  in  Figure  18.  For  the 
controlled  system,  the  set  of  feedback  gains  appearing  in 
Table  III  were  obtained  using  AMRL  Subroutine  MRIC.  The 
order  of  the  states  corresponding  to  the  feedback  gains  of 
Table  III  and  the  modal  analysis  method  is  as  follows: 

xi  = 01 
X2  = &2 
X3  = vi 
X4  = V2 
X(t)  = X3  = V3 
X6  = V4 
X7  = ui 
X8  = U2 
X9  = U3 
xicf  U4 


Ll 
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*11  = 

P01 

X12  = 

P02 

X 

II 

Pvi 

*14  = 

P V2 

^15  = 

P V3 

*16  = 

Pv4 

*17  = 

Pui 

*18  = 

Pu2 

*19  = 

PU3 

II 

0 

CM 

X 

Pu4 
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FIRST  MQIE  M0TI0N-FT  FIRST  M0IE  M0TI0N~FT 


Figure  15.  Uncontrolled  Second  Mode  Antenna  Motion 


Figure  16.  Uncontrolled  Third  Mode  Antenna  Motion 
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With  the  feedback  gains  applied  to  the  system  differ- 


ential equations  of  motion,  the  control  required  on  the 
angular  velocities  and  62  obtained.  The  resulting 
controlled  motion  of  9]^  and  02  is  depicted  in  Figure  I9. 

The  controlled  motion  of  the  first  four  modes  is  as  illus- 
trated in  Figures  20  through  23.  Figure  24  shows  the  over- 
all magnitude  of  the  antenna  motion  as  a function  of  time. 
Comparing  the  uncontrolled  to  the  controlled  motion,  that 
is.  Figure  I3  to  Figure  I9  and  Figure  18  to  Figure  24,  it 
can  be  seen  that  the  desired  null  orientation  of  the  system 
is  rapidly  achieved.  The  control  in  angular  velocity  rates 
required  to  obtain  this  orientation  is  plotted  in  Figtares 
25  and  26.  Figure  25  shows  that  the  controller  is  active 
for  only  a short  period.  Figure  26  more  clearly  illustrates 
the  magnitude  of  control  required  during  the  more  active 
stage.  This  figure  also  illustrates  that  the  controlled 
angular  velocity  rates  are  nearly  zero  after  a time  of 
one-tenth  of  a second.  If  the  actual  mechanism  imparting 
control  was  known,  these  control  angular  velocity  rates 
could  be  transformed  into  some  thrust  requirement  or  mo- 
mentum exchange  requirement. 

Recall  that  earlier  in  this  thesis  the  assumption  was 
made  that  the  results  obtained  by  the  continuous  modal  meth- 
od would  be  considered  an  exact  representation  of  the  sat- 
ellite system.  Later  in  this  thesis,  the  results  obtained 
by  applying  integral  coordinate  feedback  gains  to  the  con- 
tinuous modal  system  will  be  compared  to  the  results  just 
presented . 
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Figure  19 . Controlled  Angular  Motion  - Modal  Analysis 


Figure  20.  Controlled  First  Mode  Antenna  Motion 


Figure  22.  Controlled  Third  Mode  Antenna  Motion 


Figure  25.  Control  Requirements  - Modal  Analysis 
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Intefo-al  Coordinate  Results 

P^or  the  continuous  system,  using  integral  coordinate 
techniques,  the  uncontrolled  system  motion  is  as  depicted 
in  Figures  27  and  28.  Note  that  the  antenna  motion  is 
again  naturally  damped  out.  Additionally,  Figure  28  shows 
that  the  satellite’s  spin  axis  is  again  repositioning  it- 
self away  from  the  nominal  spin  axis  direction.  To  estab- 
lish the  control  required  to  reorient  the  satellite's 
spin  axis,  feedback  gains  were  obtained  from  Subprogram 
OPTCON.  These  feedback  gains  were  verified  by  using  the 
Subroutine  MRIC  to  determine  required  gains.  The  resulting 
feedback  gains  obtained  by  both  methods  were  identical  and 
appear  in  Table  IV. 


Table  IV. 

Integral  Coordinate  Feedback  Matrix 


Col  1 1 

Col  2 1 

col  3 , 

Col  4 1 

ow 

1 

-3.73145+02 

-3.51425+01 

-2.35395-05 

-3.25415+00 

it  ov/ 

2 

-2.35393-05 

3.03793-05 

-5.43185+02 

2.97215-06 

, Col  5 , 

Col  6 1 

Col  7 1 

Col  8 

Ro\/ 

1 

9.36715-01 

-1.66215-01 

-1.26055-01 

5.16395-01 

Row 

2 

-1.31065-01 

1 . 38205-07 

9.33475-01 

-5 . 4437-^-07 
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As  in  the  previous  techniques,  the  required  control  was 


i 

t 

! 


found  by  performing  the  matrix  multiplication  of  equation 
(55)-  For  the  integral  coordinate  method,  the  states  were 

= 9i 
= V 

= 92 

= u 

= P9i 
= 

= P92 

= Pu  (58) 

Figure  29  depicts  the  controlled  states  92_,  92,  u,  and  v 
versus  time  for  the  integral  coordinate  technique  as  ob- 
tained by  the  AFIT  Subprogram  OPPCON.  Figures  30  and  31 
illustrate  these  same  states  as  obtained  by  integrating  the 
integral  coordinate  equations  of  motion  forward  with  control 
included  in  the  system.  As  these  figures  indicate,  the 
sf  i,ellite  returns  to  the  desired  attitude  in  approximately 
one  second.  Figure  32  depicts  the  amount  of  control  re- 
quired using  integral  coordinate  techniques.  A more  de- 
scriptive illustration  of  this  control  is  given  in  Figure 
33.  By  comparing  Figure  33  bo  Figure  26,  it  can  be  noted 
that  control  requirements  for  the  integral  coordinate  and 
modal  techniques  are  nearly  identical.  The  all  negative 
eigenvalues  appearing  in  Table  V indicate  the  stability  of 
the  integral  coordinate  method  for  this  satellite  configura- 
tion. 


defined  in  the  following  orden 

^2 

^3 


x(t) 


X4 

^5 

^6 

X7 

^8 


68 


Figure  30.  Controlled  Angular  Motion  - Integral  Coordinates 
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Figure  33*  Control  Requirements  - Integral  Coordinates 


Table  V. 


Integral  Coordinate  Closed-Loop  Eigenvalues 


Real 

Imaginary 

-.47474E+01 

-.33747E+O2 

-.47474E+01 

.33747E+02 

-.56459E+01 

-.21147E+02 

-.56459E+01 

.21147E+02 

-.I8968E+03 

-.I9868E+O3 

-.I8968E+03 

.I9868E+O3 

-.27159E+03 

-.276I5E+O3 

-.27159E+03 

.276I5E+O3 

I 

r 

i Analysis 

F 

I The  results  to  this  point  have  shown  that  regardless 

I of  which  of  the  three  techniques  is  choosen  to  investigate 

i 

' the  system,  some  feedback  gains  can  be  determined  which 

return  the  modelled  satellite’s  attitude  back  to  a desired 
orientation.  Although  the  integral  coordinate  method  is 
apparently  a viable  technique  for  controlling  a satellite 
system,  further  investigation  must  be  made  to  determine  its 
effectiveness  in  providing  feedback  gains  which  could  con- 
trol an  actual  satellite  system. 

Recall  the  assumption  that  the  continuously  modelled 
modal  analysis  investigation  represents  the  physical  sat- 

; ellite  system.  Hence,  the  effectiveness  of  the  integral 

! 

I 
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coordinate  technique  can  be  determined  by  applying  integral 
coordinate  feedback  gains  to  the  continuous  modal  invest- 
igation and  seeing  if  the  satellite  does,  in  fact,  return 
to  a desired  orientation.  To  this  end,  modal  displacements 
were  related  to  integral  coordinate  displacements  by  the 
following: 


u = 


U = 


I 


h + 1 


h 


pzudz  - integral  coordinate 


2 

i = 1 


modal  analysis 


or 


Similarly, 


-h  + 1 4 


u = 


^ 2 t 

Z $:u<pzdz  = z (ml  ) S„.Ui 
h i = 1 ^ ^ . i = 1 ' ^1  ^ 


V = 


^ 2 i 

Z (ml^)^S2iVi 
i = 1 


P-  = 2u  + 2Q2  + 
RI 


P-  = ^ - 201  -f  202^ 

RI 


(59) 


Basically,  the  control  was  determined  by  first  equating  modal 
states  to  integral  coordinate  states.  After  this,  the  re- 
sulting integral  coordinate  displacements  were  multiplied 
by  the  previously  obtained  integral  coordinate  feedback  gains, 
which  yielded  the  required  control.  Numerical  integration 
of  the  now  modified  modal  system  provided  the  following 
results.  Figure  34  illustrates  that  the  satellite  does 
return  to  a desired  pointing  attitude,  as  evidenced  by 
and  02  returning  to  zero  radians.  Comparing  Figure  3^  to 
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Figure  19,  it  can  be  seen  that  the  integral  coordinate 
feedback  gains  do  control  the  system  almost  as  well  as  the 
optimal  feedback  gains  obtained  for  the  purely  modal  anal- 
ysis investigation.  Figures  35  through  38  depict  mode 
displacements,  and  again,  the  control  is  such  that  the 
antennas  return  to  a nominal  position.  Figure  39  presents 
the  magnitude  of  antenna  displacements  for  the  continuous 
modal  system  using  integral  coordinate  feedback  gains.  If 
Figure  39  is  compared  to  Figure  24,  it  can  be  noted  that 
the  magnitude  of  antenna  motion  for  the  continuous  modal 
system,  using  integral  coordinate  feedback  gains,  is  nearly 
as  controlled  as  the  pure  continuous  modal  analysis.  The 
control  required  to  return  the  satellite  to  the  desired 
orientation  is  depicted  in  Figures  40  and  41.  These  control 
results  are  comparable  to  the  results  obtained  in  Figures 
25  and  26  for  the  pure  modal  investigation.  The  above 
analysis  has  shown  that  integral  coordinate  feedback  gains 
do  control  the  satellite  system  almost  as  well  as  feedback 
gains  obtained  by  modal  techniques.  This  supports  the  sup- 
position that  the  integral  coordinate  method  is  a viable 
technique  for  estimating  optimal  control  feedback  gains  for 
an  actual  satellite  system. 

Other  example  problems  were  investigated  with  the 
results  being  basically  the  same  as  those  presented  above. 

As  long  as  the  satellite  system  was  inherently  stable, 
integral  coordinate  feedback  gains  could  be  found  and  ver- 
ified by  applying  them  to  the  continuous  modal  system. 
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For  satellite  systems  which  were  inherently  unstable,  either 
the  integral  coordinate  technique  or  the  modal  method  could 
control  the  satellite.  However,  when  integral  coordinate 
feedback  gains  were  substituted  into  the  modal  system, 
the  system  would  not  control  as  desired.  This  was  appar- 
ently due  to  the  fact  that  the  control  on  the  antenna,  which 
was  pervasive  in  nature,  was  unable  to  keep  up  with  the 
inherent  instabilities  of  the  system. 
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Figure  35.  Controlled  First  Mode  Antenna  Motion  Using 

Integral  Coordinate  Feedback  Cains  in  Modal  System 
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U2  IXSTRNCE 


Figure  36.  Controlled  Second  Kode  Antenna  Motion  Using 

Integral  Coordinate  Feedback  Gains  in  Modal  System 
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Figure  38.  Controlled  Fourth  Mode  Antenna  Motion  Using 

Integral  Coordinate  Feedback  Gains  in  Modal  Syster. 
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Figure  40.  Control  Requirements  Using  Integral  Coordinate 
Feedback  Gains  in  Modal  System 
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V . Conclusions 

3y  analyzing  the  results  to  example  1,  several  conclu- 
sions may  be  made  regarding  the  utility  of  the  three  methods 
investigated.  A discussion  of  these  conclusions  is  presented 
below. 

The  discrete  system  analysis  was  the  simplest  to  per- 
form. However,  the  results  obtained  by  this  method  are 
least  useful  to  the  satellite  designer.  The  method  was  sim- 
ple in  that  only  eight  first  order  differential  equations 
of  motion  needed  to  be  solved  to  describe  the  attitude  of 
the  satellite.  This  required  minimal  programming  effort 
and  minimal  computer  time  in  obtaining  a solution  to  the 
control  problem.  The  method  lacked  utility  in  that  the 
continuous  antennas  were  crudely  discretized  and  could  not 
be  easily  related  to  an  actual  physical  satellite.  Since 
the  model  did  not  have  the  completeness  of  the  continuous 
system  model,  a designer  would  lack  confidence  in  his  abil- 
ity to  obtain  optimal  control  feedback  gains  which  approx- 
imated those  required  for  the  true  system. 

The  purely  continuous  modal  technique  offers  a reversed 
situation  from  that  of  the  discrete  method.  The  modal  anal- 
ysis technique  offers  a means  of  obtaining  realistic  feed- 
back gains,  however,  the  method  is  rather  complex  to  easily 
employ.  Even  for  this  rather  simple  satellite,  having  only 
two  flexible  antennas,  the  method  required  the  solution  of 
the  fixed-free  beam  eigenvalue  problem  and  the  solution  of 
twenty  first  order  differential  equations.  The  computer 
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programming  effort  was  considerably  greater  than  that 
required  for  the  discrete  analysis,  and  computer  execution 
times  were  approximately  five  times  greater.  Although  time 
consuming,  this  method  adequately  models  the  satellite 
system  and  is  often  used  by  the  designer  to  obtain  feedback 
gains  and  the  associated  optimal  control  requirements. 

To  determine  control  requirements  for  a satellite  in 
the  early  stages  of  design,  a spacecraft  designer  would 
like  to  have  at  his  disposal  a technique  which  combines 
the  attributes  of  the  discrete  and  modal  investigations. 

The  results  obtained  for  the  integral  coordinate  method 
indicate  that  this  technique  could  provide  the  designer 
with  such  a design  tool.  The  method  was  generally  simple, 
requiring  the  solution  of  the  same  number  of  differential 
equations  as  the  discrete  system.  Problem  formulation  was 
further  simplified  since  the  beam  eigenvalue  problem  did 
not  have  to  be  solved.  Programming  effort  and  computer 
execution  times  were  approximately  the  same  as  the  discrete 
system  analysis.  As  shown,  the  satellite's  controlled 
motion  with  integral  coordinate  feedback  gains  applied  to 
the  continuous  modal  system  was  almost  identical  to  that 
of  the  pure  modal  method.  This  tends  to  verify  the  fact 
that  a designer  can  obtain  realistic  estimates  of  control 
requirements  by  using  the  simplified  integral  coordinate 
techniques  developed  in  this  thesis. 
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VI . Summary  and  Recommendations 


Summary 

This  th^^sis  investigated  three  methods  for  obtaining 
control  feedback  gains  of  satellites  containing  flexible 
appendages.  These  methods  included  a discrete  system 
analysis,  continuous  modal  analysis,  and  a new  integral 
coordinate  analysis.  For  each  of  the  techniques  considered, 
equations  of  motion  describing  the  satellite's  attitude 
were  developed.  Active  control  was  applied  to  the  satellite 
by  an  arbitrary  mechanism  capable  of  providing  angular 
velocity  inputs.  Control  was  required  in  the  system  to 
keep  the  attitude  of  the  satellite  at  a desired  null  orien- 
tation. A sample  problem  was  investigated  and  an  analysis 
of  the  results  for  the  three  techniques  was  performed.  The 
discrete  method  was  found  to  be  simple,  but  lacking  in 
utility  for  obtaining  realistic  control  estimates.  The 
modal  analysis  obtained  valid  control  results,  however,  the 
method  involved  rather  complex  mathematical  formulations. 

The  results  obtained  for  the  integral  coordinate  method 
indicate  that  it  is  a simple,  time  saving,  and  valid  means 
for  obtaining  estimated  control  requirements  for  satellites 
with  flexible  appendages. 

Recommendations 

Having  completed  this  investigation,  it  is  recommended 
that  the  integral  coordinate  method  be  used  for  obtaining 
control  requirements  for  the  specific  type  of  satellite 
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discussed  in  this  thesis.  This  recommendation  is  based  on 
three  reasons.  For  one,  it  is  easy  to  employ  during  early 
stages  of  design.  Secondly,  it  avoids  involved  numerical 
procedures.  And  finally,  the  effect  of  changing  various 
system  parameters  can  be  easily  assessed.  Although  the 
utility  of  the  integral  coordinate  technique  was  generally 
established  in  this  investigation,  it  is  also  recommended 
that  this  method  be  further  examined.  In  particular,  other 
satellite  configurations  or  other  spacecraft  in  general 
should  be  investigated  to  determine  the  utility  of  the 
integral  coordinate  technique  for  estimating  optimal  control 
requirements.  If  shown  effective  for  a wider  classification 


of  space  vehicles,  the  spacecraft  designer  would  realize  an 
enormous  savings  in  time  and  effort  in  obtaining  estimated 
optimal  control  requirements. 


1. 
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Appendix  A 


Kinetic  Energy  Derivation 

The  following  is  a derivation  of  the  kinetic  energj' 
expression  given  by  equation  (2). 


The  velocity  of  the  rigid  body  less  the  mass  oscilators  is 
given  by 

^B-m  " Vjj  + oo  X r (A-1) 

The  velocity  of  the  mass  oscillators  is  given  by 

Vjjj  = + (1)  X Fjjj  + ui  + vj  (A-2) 

Therefore,  the  kinetic  energy  expression  given  by  equation 
(1)  can  be  written  as 

+ i/  (V„-V„)dm  (A-3) 

B-m  m 


or 

T 


if  (Vc  + 0)  X r).(Vc  + u)  X r)dm  + if  (V^  + u)  x r^ 
B-m  m 

+ ui  + v3)«(Vjj  + (1)  X rj„  + ui  + v3)dm 


--  -T  _ __  _T_  .2.2 

= + idj  + mVjj'Vj,  + ^U)  Injd)  + m(u  + v ) 

2m(ou  X r^i)  (ul  + vj)  + terms  which  go  to  zero 

(A-4) 
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Substituting  equations  (A-5)  and  {A-6)  into  equation  (A-4) 
yields 

r = iMVc-Vc  i(Aiix^  + B.ly^  + Cilz^)  + mVc-Vc 
+ in[(a^  + "*■ 

- ZuYkl^ly  - 2aiu2x*iz  " Baviiyil^]  + 2in[uVilz  - uvQj; 

- avilx  + auilyj  + m(u^  + v^)  (A-7) 

2 

Combining  constant  terms  and  letting  A'  = A + 2ma  and 
B*  = B + 2ma^  yields 

T = i AMx^  + i 2 ni[u^  + + v^i2x^ 

-»•  u^ily^  + (u^  + v^)i22^  ■ 2a^2x  '*■  2aui2y  - 2(uv 

- uv)il2  - 2uvi]xi^y  - 2avi2yi22  - 2auA2xiizJ  Const 

(A-8) 

Equation  (A-8)  is  in  agreement  with  the  kinetic  energy 
expression  obtained  in  Ref.  4 for  the  same  satellite  system. 
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Appendix  B 

Trans format ion 

The  following  is  a transformation  from  orbitally 
referenced  angular  velocities  to  body  referenced  angles 
and  angular  rates. 


The  angular  velocity  of  the  satellite  can  he  described  by 
the  vector 

U)  = ®2®y2  (B-4) 

From  equations  (B-1),  (B-2),  and  (B-3)  it  can  be  shown  that 


e„ 

’ 

X2 

e 

y2 

_®Z2_ 

1 

0 

0 


cos  - sin 


sin  01 


cos  01 


'-3 

.^“3. 

(B-5) 


and 


®X1 

®yi 

0)  , 

. 

cos  02  sin  0]^  sin  02  cos  0i  sin  02 

0 cos  01  - sin  01 

sin  02  sin  0i  cos  02  cos  0i  cos 


®X3 

®y3 

^3- 

(B-6) 


Substituting  the  appropriate  expressions  from  equations 
(B-5)  and  (B-6)  into  equation  (B-4)  yields 

(U  = (0j^  - sin  02  ^)ex3  + (sin  0i  cos  02  63  + cos  0i  92^®y3 

• • 

+ (cos  01  cos  02  83  - sin  0i  02)ez3  (B-7) 


Hence , 


(iJ  = 


Hx 

‘V 

ilz 

10  - sin  02 

0 cos  01  sin  01  cos  02 

0 - sin  01  cos  01  cos  02 


9l 

02 

h. 


(B-8) 


94 


Appendix  C 

Derivation  of  Discrete  System  Equations  of  Motion 

The  following  is  a brief  description  of  the  derivation 
of  the  equations  of  motion  for  the  discrete  system  by  ap- 
plying equations  (13)  and  (14)  to  the  system  Hamiltonian 
and  Lagrangian.  From  equation  (13) 


^01  = = A'0i  - 2mav  - A’02kl  - 2mauA]  (C-1) 

3 0]^ 

^02  ~ 9L  = B*02  + 2mau  + B'0ii)  - C9ii2  - 2mavi2  (C-2) 

3 02 

P03  = -IL  = c93  + ca  (C-3) 

393 

Pu  = 3L  = 2mu  + 2ma92  2ma9]^il  - 2m'/il  (C-4) 

3u 

Pv  = 9L  = 2mv  - 2ma9j^  + 2ma92i2  + 2muil  (C-5) 

3v 


Equations  (C-1)  through  (C-5)  where  solved  simultaneously 
for  9j^,  9 2.  93.  u,  and  v.  From  equation  (14) 

i 

I 

= (B*  - C)i]92  + (B*  - C)a^9i  + BmsuKu  - v*])  (C-6)  | 

P02  = - A'il9i  + 2maAi(v  + m"!)  + (A’  - C)ii^92  (C-7) 

Pe^  = 0 (C-8) 

Pu  = (2mv  - 2ma9i)i2  + (2mu  + 2ma6  2)i2^  - 2Ku  - 2Du  (C-9) 

Py  = - (2mu  + 2ma92)*^  + (2mv  - 2ma9i)i2^  - 2Kv  - 2Dv(C-10) 

Substituting  in  the  relations  obtained  from  equations  (C-1) 
through  (C-5)  into  equations  (C-6)  through  (C-10),  and  after 
much  simplification,  the  expressions  of  equation  (15)  may  be 
obtained. 
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Appendix  D 


Derivation  of  Modal  Analysis  Equations  of  Motion 


The  following  is  a brief  description  of  the  derivation 
of  the  equations  of  motion  for  the  modal  analysis  method 
by  applying  equations  (37)  and  (38)*  From  equation  (37) 


Pei  = A’9i  -2  2:  ^Sz^Vi(ml^)  - A*' 


4 2 i 

- Z Sz.ui(ml'^) 

i = 1 ^ 


(D-1) 


P92  = B'e2  + 2 E Sziui(ml^)*  + B’eiil 
i = 1 


- CeiO  - 2*1  E SziVi(ml2)^ 

i = 1 


P0'j  ~ C0^  + Cd 


(D-2) 


(D-3) 


Pu-  = 2uj  + 202Szj(ml^)*  + 29idSz j(ml^)' 


- 2avj  j = 1,  — .4  (D-4) 

P^.  = 2^^  - 29iS2.(ml^)^  + 292dSz.(inl^)^ 

+ 2*luj  j = 1,  — ,4  (D-5) 

Equations  (D-1)  through  (D-5)  where  solved  simultaneously 
for  01*  02.  63.  ui,2,3,4’  ^1,2, 3, 4*  From  equation  (38) 

? 4 pi 

Pg.  = (B*  - 0002  + (B*  - C)d  01  + 20  E S2.Ui(ml^) 

i = 1 ^ 


2iT  E S2,vi(ml^)‘ 
i = 1 ^ 


(D-6) 
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(D-7) 


2 • 2 ‘ 
Peg  = -AMei  + (A*  - Oil  e2  + 2il^  l ^SziVi(ml  ) 


2 ^ 2 ^ 
+ 2^1  Z Sz.ui(ml  ) 
i = 1 ^ 


P93  = 0 


(D-8) 


Puj  = 20^3  - 2(ral2)^i^z^(ei  - 1)92)  + 2a"=Uj  - 2U)j^Uj 

[2  ^ 

^u.i  - ) (Pe,  + C9ia)  + nvj 

2 Db 


+ ^z.jtnl^  E Szj^Puj^l  j = (D-9) 

Db  i = 1 -1 


Pv  = - 2ikj  - 2(ml^)^'iSzj(92  + + 2a^Vj  - 

r ? ^ o 2 ^ 

- 2d m + ^zj(ml  ) pg  + ^zjml  z Sz.Pyi 

“pL  2 Da  Da  i = 1 ^ " 


j 


- Ouj  I j = 1.  — (D-10) 


where 


2 •*  2 

Da  = A’  ' 2ml  E Szi 

^ i = 1 

2 ^ _ 2 

Db  = B'  - 2ml  E Szi 

i = 1 


(D-11) 


Substituting  the  relations  obtained  from  equations  (D-1) 
through  (D-5)  into  equation  (D-6)  through  (D-10),  and 

after  much  simplification,  the  expressions  of  equation  (39) 

may  be  obtained. 
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Appendix  E 


T 


Derivation  of  Intepral  Coordinate  Equations  of  Motion 

The  following  is  a brief  description  of  the  derivation 
of  the  equations  of  motion  for  the  integral  coordinate 
method  by  applying  equations  (13)  and  (14)  to  the  system 
Hamiltonian  and  Lagrangian.  From  equation  (13) 


P9l 

= A'9i  - 2v  - A'92Q  - 

2ai 

(E-1) 

CM 

<1> 

= B'92  + 2u  + B'9ia  - 

C9i(l  - 2£iv 

(E-2) 

P93 

= C93  + cn 

(E-3) 

Pu 

= 2u  + 292  + 29lQ 

RI 

(E-4) 

• 

= 2v  - 29i  + 292^ 

RI 

(E-5) 

where 

RI 

^ 2 

= 1 pz^dz 

J U 

(E-6) 

Equations  (E-1)  through  (E-5)  where  solved  simultaneously 


for  9j^ 

, 92,  93,  u,  and  v.  From  equation 

(14) 

Pei 

= (B*  - C)(}92  + 

20u  + 

(B*  - 00^9 

^ - 20^v 

(E 

-7) 

P92 

= - A *091  + 2nv 

+ (A* 

- C)0^92  + 

2a^u 

P93 

= 0 

Pu 

= - 2091  + 2(0^ 

- (Ji^ 

)^  + 20^92 

• 

- 2d  u 

(E- 

10) 

RI 

pRI 

P7 

= - 2002  + 2(0^ 

- (Di^ 

) V - 2il^9i 

- 2d  7 

(E- 

11) 

RI 

pRI 
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Substituting  the  relations  in  equations  (E-1)  through  (E-5) 
into  equations  (E-7)  through  (E-11),  and  after  simplifica- 
tion, the  expressions  of  equation  (48)  may  be  obtained. 


I 

i, 

! 

I 

I 
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